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ABSTRACT 

We propose a scenario for the quasi-periodic oscillations observed in magnetar 
flares wherein a tangled component of the stellar magnetic field introduces nearly 
isotropic stress that gives the fluid core of the star an effective shear modulus. In a 
simple, illustrative model of constant density, the tangled held eliminates the problem¬ 
atic Alfven continuum that would exist in the stellar core for an organized held. For 
a tangled held energy density comparable to that inferred from the measured dipole 
helds of ~ 10 15 G in SGRs 1806-20 and 1900+14, torsional modes exist with funda¬ 
mental frequencies of about 20 Hz, and mode spacings of ~ 10 Hz. For hxed stellar 
mass and radius, the model has only one free parameter, and can account for every 
observed QPO under 160 Hz to within 3 Hz for both SGRs 1806-20 and 1900+14. The 
combined effects of stratification and crust stresses generally decrease the frequencies 
of torsional oscillations by < 10% for overtones and increase the lowest-frequency fun¬ 
damentals by up to 50%, and so the star can be treated as having constant density to a 
generally good hrst approximation. We address the issue of mode excitation by sudden 
readjustment of the stellar magnetosphere. While the total energy in excited modes is 
well within the energy budget of giant hares, the surface amplitude is < 10~ 3 of the 
stellar radius for global oscillations, and decreases strongly with mode frequency. The 
626 Hz QPO reported for SGR 1806-20 is particularly problematic to excite beyond a 
surface amplitude of 10 -6 of the stellar radius. 
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1 INTRODUCTION 


Soft-gamma repeaters (SGRs) are strongly-magnetized neutron stars with magnetic fields of B = 10 14 — 10 15 G that produce 
frequent, short-duration bursts (< 1 s) of < 10 41 ergs in hard x-ray and soft gamma-rays. SGRs occasionally produce giant 


flares that last ~ 100 s; the first giant flare to be detected occurred in SGR 0526-66 on 5 March, 1979 (Barat et al. 1979 


Mazets et al.|l979 Cline et al. 

1980 

1, releasing ~ 2 x 10 45 erg ( 

Fenimore et al.[l996 


initial peak was ~ 1 s (Hurley et al.|1999l. On December 27, 2004, SGR 1806-20 produced the largest flare yet recorded, with 


Feroci et al 


. 19991. 


The duration of the 


a total energy yield of >4 x 10 46 ergsjjln both short bursts and in giant flares, the peak luminosity is reached in under 10 
ms. Measured spin down parameters imply surface dipole fields of 6 x 10 14 G for SGR 0526-66 ( Tiengo et al.|20 09), 7 x 10 14 
G for SGR 1900+14 (Mereghetti et al. 20061, and 2 x 10 15 G for SGR 1806-20 (Nakagawa et al. 20081, establishing these 
objects as magnetars. 


* E-mail: link@physics.montana.edu 
f E-mail: anthonyvaneysden@montana.edu 
1 These energy estimates assume isotropic emission. 
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2 B. Link and C. A van Eysden 


The giant flares in SGR 1806-20 (hereafter SGR 1806) and SGR 1900+14 (hereafter SGR 1900) showed rotationally 
phase-dependent, quasi-periodic oscillations (QPOs). QPOs in SGR 1806 were detected at 18 Hz, 26 Hz, 30 Hz, 93 Hz, 150 
Hz, 626 Hz, and 1837 Hz ( Israel et al.||2005| [Watts fe Strohmayer|[2006 Strohmayer fe Watts|2006 Hnmbaxyan^et aL|2011 l. 
QPOs in the giant flare of SGR 1900 were detected at 28 Hz, 53 Hz, 84 Hz, and 155 Hz ( Strohmayer fe Watts|2005 1. Recently, 
oscillations at 57 Hz were identified in the short bursts of SGR 1806 ( Huppenkothen et al.|20T4 l, and at 93 Hz, 127 Hz, and 
possibly 260 Hz in SGR J1550-5418 ( Huppenkothen et al.|2014 l ^Overall, the data show that SGRs 1806 and 1900 support 
oscillations with spectra that begin at about 20 Hz, with a spacing of some tens of Hz below 160 Hz. The relative dearth of 
high-frequency QPOs could indicate a sparse spectrum above 160 Hz, or highly-preferential mode excitation. Given that we 
do not understand how surface oscillations affect magnetospheric emission, it could be that some modes are not seen because 
they do not create large changes in the x-ray emissivity for our viewing angle. 


Duncan (1998) predicted that the most energetic magnetar flares would excite observable seismic modes of the neutron 


star crust. The easiest modes to excite, Duncan argued, would be torsional oscillations which, for the assumed composition, 
begin at about 30 Hz. The fundamental torsional mode consists of periodic twisting of two hemispheres in opposite directions. 
Crust rigidity provides the restoring force, and the star oscillates as a torsion pendulum. 

Initial theoretical study of the observed QPOs treated the crust and core as uncoupled for simplicity ( Piro|2005| | Samuels^ 


son fe Andersson|2007 Lee|2007 Sotani et al.|2007 Watts fe Reddy|2007 Sotani et al.|2007 Steiner fe Watts|2009 ). As pointed 


out by Levin (2006), however, the strong magnetic field of a magnetar introduces an essential complication; the charged com¬ 


ponent of the fluid core, primarily protons and electrons, supports Alfven waves that couple to the highly-conductive crust 
and affect the mode spectrum. A neutron star therefore cannot be regarded as possessing crust modes that are separate from 
the core, rather, we must consider global oscillations of the crust-core system ( |Levin J200 6; Gl ampedakis et al.|2006 ). 

It has become clear that the nature of the Alfven wave spectrum of the core plays a crucial role in the global magneto¬ 
elastic problem. For smooth magnetic field geometries, such as uniform magnetization or a dipole field, the core fluid can be 
regarded as supporting an Alfven continuum of frequencies; each frequency corresponds to a natural frequency of a magnetic 
field line. Since there is a continuous distribution of line lengths through some range, the spectrum is continuous. Dipolar field 
geometries and their extensions generally give a core spectrum with continuous bands separated by gaps, and a low-frequency 
cut-off or edge, much like electronic band structure in solids. The crust, by contrast, supports only a discrete spectrum of 
torsional modes if isolated from the core. When the crust is included and coupled to the core fluid through magnetic stresses, 


the structure of the core spectrum leads to very rich dynamics as shown by Levin (2007) for a simple model with a thin crust 
and a core of constant magnetization. If energy is deposited in the crust at one of the natural frequencies of the crust, and 
this frequency lies within a portion of the core continuum, the energy is lost to the core continuum in less than 0.1 s as the 
entire core continuum is excited. The crust excitation is effectively damped through resonant absorption, a familiar process 
in MHD; 


see e.g., 


Goedbleod & Poedts (2004). 


The turning points and edges in the Alfven continuum play a special role in the dynamics; power rests at these points and 
gives QPOs. Subsequent work on purely fluid stars in general-relativistic MHD with more realistic field geometries also showed 
oscillations near the continuum edges (Sot ani et al.|2007} |2008| |Cerda-Duran etTa l. 2009; C olaiuda et al.|[2009[ |Cerda-Duran| 
et al.|2009l. van Hoven & Levin (20111 have shown that if a natural frequency of the free crust happens to fall within a gap, 


there is little loss of the crust’s oscillation energy to the core by resonant absorption; see also work by Colaiuda fe Kokkotas| 
(20111. Frequency drift, due to the effects of the continuum, generally occurs ( vanHIoven&^evin|2011|) . This basic picture 


has been further developed and refined (IGabler et al. 


jGab 

|Levin|2012| |Passamonti fe Lander|2013| [Gabier*et*aL|2013a|b|' 20141. In particular, Colaiuda & Kokkotas (2011} find that all 


2011 


Colaiuda & Kokkotas 


2011 


Gabler et al. 


2012 


van Hoven & 


observed QPOs in SGRs 1806 and 1900 can be accommodated with a particular arrangement of continuum gaps. [Colaiuda fej 


Kokkotas (2012) find in a spherical model with a magnetic field with both poloidal and toroidal components that peaks below 


about 100 Hz can persist in the power spectrum, though these peaks are far broader than observed. Gabler et al. (2013a) 


Gabler et al. (2013b), Passamonti & Lander (20131, and Passamonti & Lander (2014) find a small number of low-frequency 


QPOs that exist in the continuum gaps and at turning points for their assumed field geometries that agree qualitatively with 
observations. Superfluidity widen the gaps in the Alfven continuum ( Gabler et al.|2013b Passamonti fe Lander|[20l4 (. 

Most of the work on the global oscillation problem described above relies on tuning the locations of the gaps and turning 
points in the Alfven continuum to accommodate the observed QPOs. Here we take a rather different view than in previous 
work, and assert that the Alfven continuum that has been so problematic is unlikely to exist at all. The Alfven continuum 
results from assuming smooth, poloidal field configurations. Given the convective dynamo that is expected to operate in a 
proto-magnetar, and subsequent evolution of the initial field, we expect the field to have high-order evolving multipoles. We 
conjecture that the field is complex and tangled over length scales small compared to the stellar radius, and smooth only 
on average. High-order multipoles (the tangle), could be long lived in the core where the electrical conductivity is high, van 


^jEl-Mezeini fe Ibrahim|(2010 


by Huppenkothen et al. (2013 


reported evidence for oscillations in the short, recurring bursts of SGR 1806, but this analysis was shown 
to be flawed. 
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Torsional oscillations of a magnetar 3 


Hoven & Levin (2011) have argued that a highly-tangled field is likely to reduce or eliminate the importance of the Alfven 


continuum and to give an effective shear modulus to the core fluid; they demonstrated that the continuum is broken for a 
“box” neutron star. 

Much work on the QPO problem has included realistic neutron structure, specific magnetic field geometries, and the effects 
of general relativity. Given that the interior field geometry will likely remain unknown, we take a step back in sophistication and 
propose a minimal, illustrative model of the low-frequency magnetar QPOs (< 160 Hz). The star is essentially a self-gravitating, 
constant-density, magnetized fluid whose oscillation frequencies are determined by the anisotropic magnetic stresses from the 
large-scale organized field (primarily the dipole component) and the approximately isotropic magnetic stresses introduced 
by the tangled field; we show that density stratification and the crust typically change the oscillation frequencies by only 
~ 10% (though 30-50% for some fundamental modes), so this simple model is an adequate first approximation for developing 
a quantitative understanding of torsional modes, while elucidating the essential physics. General relativity is included only as 
a redshift factor that reduces the oscillation frequencies observed at infinity by about 20%. We show that if the energy density 
in the tangled field is comparable to that in the large-scale field, a discrete normal mode spectrum beginning at around 20 
Hz with a spacing of 10 Hz appears naturally for axial modes, generally consistent with observations of QPOs below 160 Hz, 
and that all QPOs observed in SGRs 1806 and 1900 in this frequency range can be accounted for to within 3 Hz. We do not 
address the character of high-frequency modes due to limitations to our model that we describe. To illustrate the key features 
of the model, the normal-mode analysis will be restricted entirely to axial modes. (We use the terms “axial” and “torsional” 
synonymously). 

In f|2] we derive the equations of motion for a star with a tangled field, and show how the tangled held gives the core 
fluid an effective shear modulus. In (J3| we briefly review how an Alfven continuum arises for a uniform held. In 0 we present 
analytic solutions for a star with only a tangled held; these solutions are useful for understanding the results of a more general 
held conhguration consisting of an organized held plus a tangled held, described in m in 0 we show that realistic stellar 
structure and crust stresses both have relatively small effects on most of the eigenfrequencies for torsional modes if the energy 
density in the tangled held is comparable to or larger than that in the organized held. In 0 we study the energetics of QPO 
excitation. In fj8] we discuss the implications of our work, our key results, limitations of the model, and future improvements. 
In Appendix [Aj we describe the variable separation procedure we used to obtain the solutions of 0 In Appendix [B] we use 
a simple model to show why the crust has a negligible effect. 


2 EQUATIONS OF MOTION 


The matter is subject to magnetic and material stresses. We will hnd that QPOs in SGRs 1806 and 1900 can be explained 
by average magnetic helds above 10 15 G, comparable to the upper critical held for superconductivity, so we take the star to 
be a normal conductor. The dipole held could be much smaller than the average held, and will not affect the basic picture 
provided the protons are normal; we do not consider the effects of superconducting protons in this preliminary study. 

The stress tensor for matter permeated by a held B is 


1 

4-7T 


BiBj — 



+ n (ViUj + VjUi ). 


(1) 


where /r is the shear modulus of the crust. We assume a perfect conductor, so that perturbations in the held satisfy 


SB = V x (u x B), 


( 2 ) 


where u is the displacement vector of a mass element. We specialize to shear perturbations, so that V ■ u = 0, for which eq. 
becomes 

SB = (B ■ V)it — (u ■ X7)B. (3) 


For a displacement u, the stress tensor is perturbed by 


STa = 


1 

47T 


BjBk^k'U'i BjUkVkBi 



BkBiViUk 


+ 



BkUiViB k 


+ fxViUj + transpose, 


(4) 


where repeated indices are summed and fi is the shear modulus of the matter, non-zero only in the crust, and Bi denotes a 
component of the unperturbed held. We include fi here only for completeness; we will ultimately Hnd that material stresses 
are dominated by those from the tangled held. 

We treat the magnetic held as consisting of an organized, largely dipolar contribution B 0 , plus a much more complicated 
tangled component B t \ 


B — B 0 + Bt- 


(5) 


We assume that the held is tangled for length scales smaller than It, small compared to the stellar radius. Given the un¬ 
certainties in the overall held structure, we take B 0 constant for simplicity. We denote volume averages over It as (...). We 
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4 B. Link and C. A van Eysden 


assume that different components of the tangled field are uncorrelated on average. Under this assumption, the tangled field 
can contribute only isotropic stress over length scales above It, so that 


(BiBj) = B°B° + {B^Sij, 


( 6 ) 


where (Bf) is a constant. 

To treat the tangled field’s contribution to the stress, we average the perturbed stress tensor of eq. 0- Eqs. @ and § 
to obtain 




( BjBk)VkUi — Uk(BjX7kBi) — -5ij{BkBi)Viitk + -5ijUi{B k X7 iBk) 


+ pVtiij + transpose, 


(7) 


where Ui now denotes a component of the displacement field averaged over l 3 . 

If different components of the tangled field are uncorrelated over It, one component will also be uncorrelated with the 
gradient of a different component, that is, 


(BiVkBj) = ((B° + Bt)V k (B? + B))) = =0 i^j. 


( 8 ) 


Since the tangled field varies over length scales smaller than It, a component of the tangled field will also be uncorrelated with 
the gradient of the same component, so that 

{B\V k B\) = 0, (9) 

which applies component by component, and therefore also in summation, as given above. Using eqs. 0, 0- 0, and V--u = 0, 
eq. 0 becomes 

(STij) = -j— ( Bj B k S7 k Ui + B t B k X7 k Uj — 5ijB k Bi V mAj + {Bt) + (V tUj + V^tti). (10) 

The tangled field gives the fluid an effective shear modulus of (B?)/4n, and enhances the rigidity of the solid. 

Upon comparing our mode calculations with data, we will find that the total energy in shear waves in the core greatly 
dominates that in the crust. We henceforth ignore crust rigidity, and justify this approximation in S0and Appendix [TT| 

We will be interested in modes with wavelengths greater than It, for which the equation of motion is 

d 2 Ut 


Pd- 


dt 2 


= Vi(5Tjj), 


(ID 


where pd is the dynamical mass density of matter that is frozen to the magnetic field. If the the protons are normal, as we 
have assumed, there will no entrainment between the protons and neutrons; entrainment of the neutrons is in any case a small 
effect if the both the protons and neutrons are superfluid (Charnel & Haensel 2006). We take pd to be the proton density x p p, 
where p is the mass density, and x p ~ 0.1 is the average proton mass fraction in the core. 

We neglect coupling of the stellar surface to the magnetosphere, and treat the surface as a free boundary with zero 
traction, thus ignoring momentum flow into the magnetosphere. Under this assumption, the traction at the surface vanishes: 


ri{STij) = 0, 


( 12 ) 


where f is the unit vector normal to the stellar surface. 

Given the uncertainties in the field geometry, we henceforth consider a uniform star of density p = 3M/4 ttR 3 for 
illustration, where M and R are the stellar mass and radius, and take B 0 = zB 0 where B 0 is constant. Eqs. (101 and (111 
give 

S 

° dz 2 

where (? a = B 2 /(4npd) and c 2 = (B 2 )/(4npd)', c 0 is the speed of Alfven waves supported by the organized field, and Ct 
is the speed of transverse waves supported by the isotropic stress of the tangled field. We show in 0 that realistic stellar 
structure changes the eigenfrequencies of torsional modes by typically ~ 10% (but by up to ~ 50% for the lowest-frequency 
fundamentals), so that a constant-density model is a good first approximation. 


2 d~U 2t -7 duz 2vy2 . 2 n 

c„ -rvr — c„V — - 1- Ct V u + ui u = 0, 

dz 


(13) 


3 THE ALFVEN CONTINUUM 


For finite B a and B t = 0, there exists a continuum of axial modes (Levin p007| , given by eq. © 
du z 


2 d'u 2v7 

7 o C oV , 

dz 2 dz 


+ UJ U = 0, 


(14) 


In cylindrical coordinates (s,<j>,z), axial modes are given by u = For a constant field, and no crust, field lines have 

a continuous range of lengths between zero and 2 R, determined by the cylindrical radius s. Within the approximation of 
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Torsional oscillations of a magnetar 5 


ideal MHD, fluid elements at different cylindrical radii cannot exchange momentum. At a given s, the length of a field line is 
2V R 2 — s 2 . The solutions have even parity (cos fez) and odd parity (sin fez). The requirement that the traction vanish at the 
stellar surface gives the spectrum 

Co 


n 

Vn = - 


2 VR 2 ~ s 2 
(2ra + l) 


even parity 


VR 2 ~ . 


odd parity. 


(15) 


(16) 


where n is an integer, beginning at zero for the odd-parity modes, and v = u}/2ir. Because s is a continuous variable, for 
every n there is a continuous spectrum of modes for this simple magnetic geometry. An infinite sequence of continua begins at 
~ 7(2n + 1) B° 5 Hz for odd-parity modes and ~ 7(2n) B° 5 Hz for the even-parity modes, where B° 5 = B 0 / (10 15 G). The full 
spectrum begins at ~ 7 B° 5 Hz; there is also a zero-frequency mode corresponding to rigid-body rotation. The same conclusion 
holds for more general axisymmetric field geometries, though certain geometries give gaps in the continuum. 


4 ISOTROPIC MODEL 


We now turn to the opposite extreme of a tangled field that dominates the stresses, taking B 0 = 0, and solving the resulting 
isotropic problem. This problem provides useful insight into the mode structure of the more general problem with non-zero 
B 0 and (B 2 ). For this case, eq. (131 becomes 

2 t -7 2 . 2 _ 

Ct V U + LU u = 0 , 

Subject to the restriction V ■ u = 0, the solutions for spheroidal modes (u r = 0), can be separated as 


u.$ = w{r)-^Yim{9,4>)e l 


ue = —w(r) — 


d 




sin 6 dcf> 

The radial function w(r) satisfies Bessel’s equation: 
d 2 2 d l{l 1) 2^\ / \ r, 

^ + -y r ~^^ + k NM = o, 


(17) 

(18) 

(19) 

( 20 ) 


where fc = uj/ct- 

The solutions that are bounded at r = 0 are the spherical Bessel functions ji(kr). Zero traction at the stellar surface 
gives 

'4h_3i 

. dr r 


= 0 


( 21 ) 


For each value of l, eq. (211 has solutions xi, n = ki,„R, where n = 0,1,2..., the overtone number, gives the number of 


nodes in ji(r). The eigenfrequencies are 


&l,n — [ 1 


Rs(M) 
R 


(Bt)R 
3 x p M 


Xl,-, 


( 22 ) 


where a redshift factor z = (1 — R S (M)/R ) 1 ^ 2 has been introduced; R a is the Schwarzchild radius. In terms of fiducial values 


Vl ’ n (HZ) ~ fcr" “ 4 ' 3 ( 0.77) (10 km) ( 1. 


M 


4M(7) 


- 1/2 


(S) 


- 1/2 


{si) 


2 \ 1/2 


10 15 G 


Xln Hz. 


(23) 


In Table [lj we give some of the eigenfrequencies for these fiducial values. We note that for {B 2 } 1 ^ 2 = 10 15 G, this simple 
model gives fundamental frequencies below 20 Hz, with a spacing of about 15 Hz for the overtones. 

For l = 1, eq. (211 has a solution w(r ) = r for k — 0. This solution corresponds to rigid-body rotation and we label it 
n = 0. This solution is of no physical significance to the mode problem we are addressing, but we include it in Tables 1-3 for 
completeness. 

Frequencies above 160 Hz correspond to k Ut iR> 30, or wavelengths <0.2 R. At these high wavenumbers, the wavelength 
of the mode could be comparable to or smaller than the length scale It over which the field can be considered tangled, in 
which case the averaging procedure of n would break down. In this case, the mode frequencies will be determined, at least 
in part, by the detailed (and unknown) field geometry. Henceforth, we restrict the analysis to v < 160 Hz, the range in which 
the low-frequency QPOs lie. 
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6 B. Link and C. A van Eysden 


1 

71 — 0 

n = 1 

n = 2 

n = 3 

72 — 4 

72 — 5 

72 = 6 

1 

0 

24 

39 

52 

66 

79 

93 

2 

li 

30 

45 

58 

72 

85 

99 

3 

16 

36 

50 

64 

78 

92 

105 

4 

22 

41 

56 

70 

84 

98 

111 

5 

27 

46 

61 

76 

90 

103 

117 

6 

31 

52 

67 

81 

95 

109 

123 


Table 1. Eigenfrequencies in Hz for the fiducial values of eq. (]23[). 


5 ANISOTROPIC MODEL 


We now turn to the general problem of non-zero B 0 and (R 2 ), and show that even a small amount of stress from the tangled 
field breaks the Aflven continuum very effectively. The normal modes of the system are similar to those found in the isotropic 
problem for (Bt) — B 2 . 

The modes are given by eq. (131. This vector equation, subject to boundary conditions on the surface of a sphere, is 
difficult to solve in general. For illustration, we specialize to axial modes, ue = u r = 0, giving in cylindrical coordinates 

2 2 t -7 2 2U<t> .2 n , nA \ 

c °-^ + c t v u^, - Ct-g^ + uo =0. (24) 


The third term results from derivatives of (j> in the original vector equation (131. Defining b 2 = (Bf)/R 2 , the ratio of the 
energy density in the tangled field to that in the organized field, the above equation becomes 

<9 2 W0 


,21 d ( du^\ / ; 


_ u 2V ± 

dz 2 4 s 2 


+ bl U ( j > = 0 . 


The zero-traction boundary condition (eq. 12 ) at the stellar surface is 


( 


du</, 


ds 




4^=o. 

oz 


#■) + (! + & t) ' 

Equation ( |25[ ) is solved in the domain 2 > 0, and so at 2 = 0 we require 

Ufj) — 0 , 


(25) 

(26) 

(27) 


for modes with odd parity about 2 = 0, and 


dz 


= 0, 


(28) 


for modes with even parity. The system is separable with a coordinate transformation; the details are given in Appendix 
|A| We solve eqs. (A8l, (A9l, and (A101 numerically to obtain the normal mode frequencies. The solutions are given by two 
quantum numbers: k and the overtone number n. k maps smoothly to 1(1 + 1) in the limit c 2 —> 0, so we use l and n for 
convenience in labeling the modes. We refer to n = 0 for a given l as the fundamental for that value of l. 

The mode structure is shown in Figure[l]for constant total magnetic energy, with (Bg + (R 2 )) 1 / 2 fixed at 10 15 G. Modes 
for which u$ has even parity (odd l) and odd parity (even l ) about 2 = 0 have been plotted separately. For large b 2 the 
isotropic solution presented in Table [l] is recovered. As b 2 is reduced, the modes become more closely spaced, approaching a 
continuum as b 2 —¥ 0. In the limit bt —> 0 the sequence of continua begins at ~ 7(2 n + 1) Hz for odd-parity modes (even l) 
and ~ 14n Hz for even-parity modes (odd l), in agreement with the continuum sequences described by eqs. (151 and (161. 
The n = 0 odd modes approach zero as 6 2 —¥ 0 and scale as 

- 1/2 


v = 6.70 1 bt 


fe)( 


R 


10 km 


1/2 


M 


1.4 M q 


(S) 


M' 1/2 ( 


B 0 


VlO 15 G, 


Hz. 


(29) 


These modes approach rigid-body rotation solutions of the continuum coupled together by the tangled field in the bt —> 0 
limit. Examination of the eigenmodes as bt —> 0 also shows that the oscillation amplitude becomes sharply peaked at a specific 
value of the cylindrical radius s and vanishes everywhere else, in agreement with the continuum solution in (|3] 

The splitting of the Alfven continuum is shown in a different way in Figure [2] where realistic predictions are made for 
the eigenfrequencies in SGR 1806 (left) and 1900 (right). In contrast with Figure [l] where the total magnetic field was held 
constant (( Bg + (R 2 )) 1 ' 42 = 10 15 G), we now fix the ordered field B a in each magnetar to its observed dipole spin-down 
value: 7 x 10 14 G for SGR 1806 and 2 x 10 15 G for SGR 1900. We take the fiducial values of M = 1.4 Mg, R = 10 in 
each magnetar. Eigenfrequencies for all n and l for the range of frequencies shown are plotted. For bt = 0, there exists an 
Alfven continuum that begins at ~ 7 R° 5 Hz, with a zero-frequency solution corresponding to rigid-body rotation. As 6 2 is 
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Torsional oscillations of a magnetar 7 




Figure 1. Splitting of the Alfven continuum as a tangled field is added to a smooth field. We have fixed (Bg + (Bj)) 1 ' 2 = Bo(l + bj ) 1 / 2 
to 10 15 G to illustrate the smooth transition from the continuum to the is otro pic tangle. For l odd and n = 0 (upper left panel), 
low-frequency modes exist that go to zero frequency for bj -4 0; see equation (|29|). 
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Figure 2. Splitting of the Alfven continuum (the column on the far left of each figure) as a tangled field is added to a smooth field in 
SRG 1806 (left) and 1900 (right). The ordered field is fixed at the dipole spin-down value: 2 X 10 14 for SGR 1806 and 2 X 10 15 for SGR 
1900. Red-dashed lines correspond to observed QPO frequencies. Solid-blue squares represent the value of bj that best fits the data in 
each object, though note that only three observed QPOs are shown on the scales of these plots. The lowest sequence of points approach 
the k = 0 rigid-body mode (not shown here), according to eq. (|29[) . 


increased, the continuum splits into discrete normal modes. As bj is further increased and the modes spread out further, the 
higher-frequency modes move off the diagram. For large bj the eigenfrequencies become large due to the large tangled field. 

The observed QPO frequencies in SGRs 1806 and 1900 are shown in Figure [2] as horizontal red dashed lines. For given 
M and R, the model has only bj as a free parameter, and gives quantitative predictions. Assuming the lowest observed QPO 
corresponds to the lowest eigenfrequency in the spectrum fixes the value of bj in each object, giving 5.0 and 0.87 for SGR 
1806s and 1900 respectively. Eigenfrequencies for these bj are plotted as solid blue squares. 

In Tables [ 2 ] and [ 3 J we give frequencies for SGR 1806 and 1900 at bj = 5.0 and 0.87 respectively. Numbers in boldface 
denote normal mode frequencies that are within 3 Hz of an observed QPO in both SGRs 1806 and 1900, and represent 
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1 

n — 0 

n = 1 

n = 2 

n = 3 

n = 4 

n = 5 n = 6 

1 

0 

39 

62 

85 

108 

130 153 

2 

18 

49 

73 

96 

119 

142 

3 

27 

58 

81 

104 

126 

148 

4 

36 

67 

91 

114 

137 

160 

5 

44 

76 

100 

123 

145 


6 

52 

84 

109 

133 

155 


7 

59 

93 

118 

142 




Table 2. Eigenfrequencies in Hz for SGR 1806 with B 0 = 7 X 10 14 and if = 5.0. Numbers in boldface are within 3 Hz of an observed 
QPO, and represent plausible mode identifications. We do not list frequencies above 160 Hz. 


l 

n = 0 

n = 1 

n = 2 

n = 3 

1 

0 

52 

89 

126 

2 

28 

71 

108 

144 

3 

38 

75 

109 

145 

4 

50 

91 

127 


5 

60 

98 

130 


6 

70 

112 

147 


7 

81 

121 

152 


8 

91 

134 



9 

101 

144 



10 

111 

155 




Table 3. Eigenfrequencies in Hz for SGR 1900 with B 0 = 2 X 10 15 and bf = 0.87. Numbers in boldface are within 3 Hz of an observed 
QPO, and represent plausible mode identifications. We do not list frequencies above 160 Hz. 


plausible mode identifications. Within these tolerances, every observed QPO frequency below 160 Hz in SGRs 1806 and 1900 
is accounted for. For SGR 1806, the model does not give the distinct 26 Hz and 30 Hz that have been reported, but we note 
that these features are quite broad (Watts & Stroh mayer|[20 06). The reported frequencies have uncertainties, and drift with 
time, and analysis has not yet been done to determine if the data favor two frequencies around 26-30 Hz, or only one. In 
our model, both are accommodated by the 27 Hz frequency for this example. If distinct QPOs at 26 and 30 Hz could be 
ascertained, this would challenge our model. 

Examples of eigenmodes are plotted in Figure [3] The left panel shows the 93 Hz mode of SGR 1806, with the attribution 
l = 7, n — 1, while the right panel shows the 155 Hz of SGR 1806, with the attribution l = 10, n = 1. The left panel has 
$ = 5.0, which is close to the isotropic limit, and shows a large degree of spherical structure symmetry. In the plot on the 
right, fe? = 0.87, and the mode shows mostly cylindrical structure, a characteristic feature of the continuum. 



Figure 3. Example eigenmodes for SGR 1806s (left) and 1900 (see text for details). Blue denotes motion out of the page for a given 
phase, red denotes into the page, and white denotes zero amplitude. The left panel has = 5 and exhibits the spherical structure of the 
isotropic limit, while the right has if = 0.86 and exhibits the cylindrical structure of the continuum. 
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1 

n = 0 

n = 1 

n = 2 

n = 3 

n = 4 

n — 5 

n = 6 

1 

0 

-15% 

-15% 

-16% 

-15% 

-16% 

-16% 

2 

+48% 

-8.9% 

-12% 

-14% 

-15% 

-15% 

-15% 

3 

+41% 

-4.5% 

-10% 

-12% 

-13% 

-14% 

-14% 

4 

+35% 

-1.9% 

-8.0% 

-11% 

-12% 

-13% 

-14% 

5 

+31% 

-0.47% 

-6.5% 

-9.4% 

-11% 

-12% 

-13% 

6 

+27% 

0.33 % 

-5.6% 

-8.5% 

-10% 

-11% 

-12% 


Table 4. Percentage frequency change of modes in a stratified star with a crust compared with a constant density star with the same 
mass and radius. We use a 1.4 Mq star with radius 13.1 km, taking x p = 0.1, and magnetic field ( B 2 ) = 10 15 G 


Though our simple, constant density model is intended mostly for illustration, it accurately reproduces all observed 
QPOs in SGRs 1806 and 1900 though the spectrum is rather dense with a spacing of about 10 Hz. The main goal has been 
to show in a quantitative manner how the Alfven continuum is broken. As we show in the effects of realistic stellar 
structure are generally small (except for low-Z fundamentals), and so this constant-density model is sufficiently accurate for 
semi-quantitative study of torsional oscillations. Many of the predicted frequencies have not been observed, and we discuss 
this further below. 

All examples shown here are for R = 10 km. If R is reduced, it is necessary to decrease b 2 to obtain similar fits to the 
data. The spectrum becomes generally less dense. Larger R has the opposite effect. 


6 EFFECTS OF STRATIFICATION AND THE CRUST 


We now show that density stratification and crust rigidity gives corrections to the eigenfrequencies of the constant-density 
model considered thus far of typically ~ 10%, but up to 50% for low-Z fundamentals. Because our best fits to the data for 
the constant-density model give b 2 = 5 for SGR 1806 and b? ~ 1 for SGR 1900, we ignore the organized field, and study the 
isotropic problem corresponding to a strong magnetic tangle. 

We construct a relativistic star using the analytic representation of the Brussels-Montreal equation of state derived by 


Potekliin et al. (2013). For a 1.4 Mg neutron star, we obtain a radius of 13.1 km. We take x p = 0.1, and and fix the magnetic 


field at (B 2 ) = 10 15 G. The shear modulus is from Strohmayer et al. (1991). For this spherical problem, the variable separation 


proceeds as in SJ4] We solve numerically for the radial eigenfunctions and eigenfrequencies. As a baseline for comparison, Table 
[I] can be scaled according to eq. (231 to obtain the eigenfrequencies for a constant-density star of radius 13.1 km. We assume 
that the dynamical density in the crust is equal to the total mass density, so that the entrainment of superfluid neutrons in 


the crust Charnel (2005 2012) is effectively perfect. Entrainment, or lack of entrainment, has little effect on our results. 


To explore the individual effects of stratification and the crust, we first examined a stratified star with no crust. We 
find that stratification lowers the frequencies for all modes with the exception of the n = 0 modes, which are raised; typical 
changes are < 20% for n = 0 and -15% for n > 0. This general trend was also found by Passamonti & Lander] |2014|. The 
reason for the frequency decrease (except for n = 0) can be seen by noting that the frequency scales as R 1 ^ 2 in the case of 


uniform density; see eq. (231. The stratified star has more mass concentrated in the central regions of the star, making it 


behave effectively as a smaller star compared to a constant-density star of the same radius. 

Next, we considered a stratified star with a crust. The addition of a crust increases all frequencies by ~ 5 — 20% for n = 0 
and 1-5% for n > 0 compared with the stratified case. The crust increases the frequency because it increases the average shear 
modulus of the star. The effect is small because the shear modulus of the crust exceeds (B 2 ) /4 -7T in only a small fraction of 
the stellar volume. The effects of crust rigidity are therefore small compared to the rigidity introduced by the tangled field 
throughout the star, and the Z = 1 modes are unchanged to two significant figures. A simple analysis is given in Appendix [B] 
that further quantifies this effect. By contrast, the rigidity of the crust is essential an essential ingredient if the field is not 
significantly tangled, as considered in previous work ( e.g ., van Hoven fe Levin|2011 Colaiuda fe Kokkotas|20TT Gabler et al. 


|2011[ |2012| |van Hoven fe Levin|2012| |Gabler et al.|2013a|b 2oT^~ 

On the balance, stratification has a larger effect than does the crust. The combined effects of stratification and a crust 
is presented in Table [4] where the percentage change in frequencies compared with a constant density star of 13.1 km radius 
are given. For the n = 0 modes, the combined effects of stratification and a crust increase frequencies by up to 50%, because 
both effects add. For n > 1 the frequencies drop as a result of stratification, but are enhanced slightly by the crust giving a 
net change of up to 16%. Changes become more pronounced for increasing n, but less pronounced for increasing Z. 

We conclude that for star in which the tangled field accounts for much of the magnetic energy density (b 2 >1), that 
both density stratification and the crust rigidity have a relatively small effect on the frequencies of torsional normal modes. 
The constant-density model is an accurate first approximation at the 10% level for most overtones, and better than 50% for 
fundamental modes. 
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7 ENERGETICS 


The only quantitative discussion to date of the energetics of QPO excitation has been given by Levin & van Hoven (2011). We 


use the results of their analysis to calculate the mode amplitude at the stellar surface in the global oscillation interpretation 
of QPOs, and the total energy in the low-frequency QPOs. 

The natural candidate for the energy source that excites stellar modes is the flare itself. For giant flares, the energy 
cannot be released deep inside the star, as the impedance mismatch between the stellar interior and the magnetosphere is so 
great that the energy in Alfven waves would take seconds or longer to reach the magnetosphere due to multiple reflections, 
in conflict with observed rise times of < 10 ms (Link||20141. This constraint favors an external origin for flares, as suggested 
by Lyutikov (2003, 2006), Komissarov et al. ( 20(*) l 7| , and Gill fe Heyl (20101. Forcing of the magnetosphere through episodic 
release of internal magnetic energy could produce a magnetospheric explosion when an MHD instability is triggered. 


Levin & van Hoven (2011) considered the excitation problem under the assumption that the magnetosphere undergoes 
a sudden change (over several relativistic Alfven wave crossing times in the magnetosphere, ~ 10 -4 s) from one equilibrium 
configuration to another, exciting torsional oscillations in the star. Suppose the shear stress at the surface changes by fB 2 /4n , 
where B is the characteristic field strength in the inner magnetosphere, and / < 1. The average displacement Ui at the stellar 
surface, assuming excitation of global mode i with frequency is of order (Levin & van Hoven|2011 


eq. 11) 


Ui 

~R 


1 fB 2 R 3 fE„ 


4n Men 2 2tt MR 2 Lof 


= 1CT 


fEu 


2 x 10 47 erg 


J \10km) \1 AMqJ V30 Hz/ 


(30) 


follows from the Fourier transform of 


where E mag = (4nR 3 /3)(B 2 /8 tt) ~ 2 x 10 4 'Bi 5 erg is the characteristic magnetic energy available to be released; /F mag is the 
total energy released, and should be comparable to the flare energy. The scaling as oj~ 
a step function. 

The energy in a mode is of order ( |Levin fe van H oven|201 l[ eq. 12) 
f 2 B 4 R 4 


E; 


327r 2 Mcn 2 ' 


(31) 


Most of the energy is in low-frequency modes. We have found a fairly dense spectrum with a characteristic mode spacing of 
Av ~ 10 Hz. The sum of the total energy in all modes can be approximated with an integral 


E ~ 


du 


Eiy) 

Av 


10 44 f 


B 


(lO 15 G 


V f^r 1 f-^-v 

) V10 km/ \1.4MqJ \20 Hz/ VlOHz/ 


erg, 


(32) 


where vo is the lowest-frequency normal mode in the spectrum. We see that the conversion of magnetospheric energy to 
mechanical motion in the star an inefficient process; about 0.1% of the flare energy goes into normal modes. The total energy 
in the modes in our model is well within the energy budget of a giant flare. 

10 8 than if _E maE is liberated. These 


For a small flare, with /E mag = lCr a erg, the amplitude from eq. 30 is smaller by 


estimates hold not just for our model, but generally for any global-oscillation model of QPOs. 


8 DISCUSSION AND CONCLUSIONS 

We propose a new physical scenario for magnetar oscillations wherein the Alfven continuum supported by the organized field 
is broken by nearly isotropic stresses arising from components of the field that form a complex tangle over scales much smaller 
than the stellar radius. In a simple model of constant density, an organized field of strength ~ 10 15 G, and a tangled field of 
comparable strength, a fundamental frequency appears near 20 Hz for 1 = 1 and l = 2, with a with a mode spacing of ~ 10 Hz 
(we use l from the solutions of the isotropic problem to label the modes of the anisotropic problem). The model is consistent 
with QPOs observed in SGRs 1806 and 1900 below about 160 Hz. Magnetic stresses dominate material stresses almost 
everywhere in the crust, so crust elasticity is unimportant when a strong, tangled field is present. We find that the combined 
effects of realistic stellar structure and crust stress give a modest change in the eigenfrequencies of ~ 10% for most modes 
(but up to 50% for low-Z fundamentals) so our simple model with constant density is a good first approximation to studying 
tangled fields that elucidates the essential physics. The strength of the organized field, which should be approximately dipolar, 
is determined observationally. Once stellar mass and radius are fixed, the model has only one important free parameter: the 
ratio of the energy densities in tangled and organized fields. 

Our results are insensitive to whether or not the protons are superconducting. The upper critical field for type II 
superconductivity is < 10 16 G, comparable to average fields we find to be consistent with the QPOs observed in SGRs 
1806 and 1900; the magnetic stress is nearly the same if the protons are superconducting, though we have assumed for 
definiteness that they are normal. Also, if the protons are superconducting, the neutrons are only slightly entrained by the 
protons throughout most of the core, and the dynamical mass density is nearly equal to the proton mass density in this case 
as well. 
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Torsional oscillations of a magnetar 11 


While our constant-density model is mostly illustrative, it is nevertheless sufficiently accurate to study torsional oscillations 
semi-quantitatively within our approximations; the model is able to account for every QPO observed in SGRs 1806 and 1900 
to within 3 Hz. No other published model has given such quantitative agreement, though we acknowledge that the success 
of the model is due in large part to the rather dense spectrum we predict with a frequency spacing of about 10 Hz. We find 
the data are best explained if there is approximate equipartition of the energy in the smooth and tangled fields. Such modes 
have high enough wavenumbers to probe the field structure over scales below 0.2 R, at which point the averaging procedure 
we introduced to treat the tangled field might become inappropriate. 

Why have most of the predicted frequencies not been observed? One possibility is that the modes have been excited, but 
are not visible. Given that we do not understand how surface oscillations affect magnetospheric emission, it could be that 
some modes are not seen because they do not produce sufficiently large changes in the x-ray emissivity to be visible from our 
viewing angle. Another possibility is that the instability that drives the flare creates highly-preferential excitation of axial 
modes. We note that preferential excitation of low-frequency free oscillations occurs in the Earth (Rhie & Ro manowicz|200~4 1, 
the so-called “hum”. The question of how preferential excitation of stellar modes might occur is an interesting question. 

While our focus has been on SGRs 1806 and 1900, which have dipole fields of ~ 10 15 G, QPOs have been reported for 
SGR J1550-5418 (Huppenkothen et al.|2014) in association with burst storms. SGR J1550-5418 has an inferred dipole field 
of 3.2 x 10 14 G ( Carnilo et al.|2007 |. In this object, as in SGR 1900, the observed QPOs can be accounted for if most of the 
magnetic energy is in the tangled field, so that the average magnetic field is ~ 10 15 G. 

A crucial ingredient in the interpretation of QPOs as stellar oscillations is to understand how crust movement can produce 
the large observed modulations of the x-ray emission by 10-20%. Timokhin et al. (2008) propose that twisting of the crust, 
associated with a stellar mode, modulates the charge density in the magnetosphere, creating variations in the optical depth 
for resonant Compton scattering of the hard x-rays that accompany the flare. In that model, the shear amplitude at the 
stellar surface must be as large as 1% of the stellar radius. D’Angelo & Watts (20121 accounted for geometrical effects of 
the beamed emission, and find that the amplitude of the QPO emission is increased by a factor of typically several over the 


estimate of Timokhin et al. 2008) for a given amplitude of the surface displacement. Based on the analysis of Levin & van 


Hoven (2011) of mode excitation by a sudden change in magnetospheric equilibrium, we find that the predicted, low-frequency 

44 erg in mechanical energy. This number is well within the energy budget of giant 
3 1? (30 Hz/V) 2 , smaller than required by Timokhin et al. (20081, especially for 


QPO spectrum could contain up to ~ 10 3 
flares but the surface amplitude is only 


10 " 


the QPOs near 160 Hz. This small amplitude could spell trouble for the global-oscillation interpretation of QPOs, unless a 

-.39 


more efficient mechanism than that of Timokhin et al. (20081 can be identified. Explaining QPOs in small bursts (~ 10 3 
erg) appears particularly difficult with a global oscillation model. The mode amplitude scales as the flare energy (eq. 301, 
giving a surface amplitude that is ~ 10 s times smaller than for a giant flare. The 626 Hz QPO reported in SGR 1806 also 
seems problematic to excite in a giant flare, with a surface amplitude of ~ 10 _6 i? in the most optimistic case. These numbers 
show that energetics should be given serious consideration in QPO models, and that other excitation mechanisms than that 
of Levin & van Hoven ( 2011) should be considered. 

The simple model presented in this paper can be improved by including a crust and introducing realistic stellar structure. 
We have shown that the effects typically ~ 10 for most modes for a star in which the magnetic energy density is dominated 
by that in the tangled field; the effects of stellar structure and the crust will become more important as the effective shear 
modulus from the tangled field is reduced, and we have not quantified these effects. Other modes than axial modes, such as 
polar modes, should be considered as well. 

An explanation for why the observed frequencies are quasi-periodic is well beyond the scope of our model. Quasi-periodicity 
could arise from magnetospheric effects, dynamics in the star not included here, or both. 
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APPENDIX A: VARIABLE SEPARATION 


Eqs. (251 and ( |26| ) can be separated and solved by transforming to an oblate spheroidal coordinate system ( u , v) defined by 

(Al) 
(A2) 


s = R 
z = Ruv 


(1 + b 2 u 2 ) (1 — v 2 ) 

i+6 2 
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Curves of constant u are ellipses, and curves of constant v are hyperbolae. For it = 1, the coordinate gives a sphere of radius 
R. In the limit 6? —> oo, spherical coordinates are recovered with u = r/R and v = cos (9. 

In these coordinates, eq. (25 I becomes 
=>2„ 


/, . t 2 2 \ b) ^(f> . ri 12 btUc/} ,2 / 2 \ & 

( 1 + fe ^)^F+ + 


0 ,2 9 u$ 0 t u,f, _ 2 / 2 2 . 2\ n 

-2b t v-g^ - - _ + UJ [b t u + v)u</ > = 0 


b t u 0 


where 

CJ = - 

Ct 


b ? 


1 + 6 ?’ 


and c? = {Bt)/(4np d ). 

The boundary condition eq. (261 at u = 1 becomes 
(1 + fe?)^-^ = 0 , 
while at a = 0 we require 

U(ji — 0, 

for odd parity modes or 

(l + b 2 t u 2 ) v ^+b 2 t (l-v 2 )u^=0, 

for even parity modes. 

We seek a separable solution of the form 

u^(u,v) = U(u)V{v). 

Eq. (|A3|) becomes 


- Kb 2 U + u) 2 b 2 u 2 U = 0 


h , ,2 2X d 2 U 2 dU 

^1 bfU j ——— H - QibfU —— H - 


du 2 


du 1 + 6?m 2 


,2 2\ d 2 V 2 dV b 2 v ,2t. _2 2 T/ n 

b t U - v ) -—it - 2b t v— --- ? + nb t V +uj v V = 0, 

' ' av A dv 1 — v 2 


where k is the separation constant. The boundary condition eq. (A4l becomes 

d + ^f-6?!/^, 


at u = 1. The boundary conditions (A5l and (A6) are satisfied if we impose 

U{ 0) = 0 and E(0) = 0, 

for odd parity modes or 

U’( 0) = 0 and I/'(0) = 0, 


(A3) 

(A4) 

(A5) 

(A6) 

(A7) 

(AS) 

(A9) 

(A10) 

(All) 

(A12) 


for even parity modes. 

In the limit 6? —> oo, eq. (A8) reduces to the spherical Bessel equation with solution U(u) = j(um), while (A9l reduces 
to an associated Legendre equation with solution V(v) = \/l — v 2 dP°/dv = —Pl(v) and k = 1(1 + 1), which recovers the 
solutions of the isotropic model of [j4] 

We solve eqs. ( |A8[ ), ( |A9[ ), ( |A10[ ) numerically for the normal mode frequencies. 

For arbitrary 6? the modes are identified by the number of nodes; for the n th mode the function U(u) has n nodes on the 
domain 0 < u < 1, while for the I th mode the function V(v) has l — 1 nodes on —1 < v < 1. 


APPENDIX B: EFFECTS OF THE CRUST 

The comparisons to observed QPOs in fJHJsuggest 6? ~ 1, so that (Bf)/4n ~ B 2 /4n ~ 8 x 10 28 erg cm -3 . The shear modulus 
crust exceeds this value only in the densest regions of the crust. The rigidity of the crust will increase the normal-mode 
frequencies somewhat with respect to what we have found by neglecting the crust. Here we show that crust rigidity has small 
or negligible effects on our results. 

To estimate the effect, we consider a two-component isotropic model with only a tangled held. The core liquid, which we 
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take to be homogeneous, has an effective shear modulus fj, co re = (B 2 ) /47 r and a shear-wave speed Ct . The crust, which we also 
take to be homogeneous, has shear modulus /r cr ust and a shear-wave speed e M . The crust has an inner radius R c and outer 
radius R. 

In the core, the solution to the mode problem is u c ore = ji(kr); the mode frequency is ai = Ctk. The solution in the crust 
is itcrust = aji(k'r) + bni(k'r), where ni are spherical Neumann functions, a and b are constants, and u = c^k'. The boundary 
conditions are continuity in value and traction at r = R c , and vanishing traction at r = R\ 


Ucore ( Rc ) — Ucrust (Rc ) , 

du C c 


/^C 


U c , 


. dr 


l^c rust 


C^crust derust 


dr 


— /^crust 


= 0 . 


du c 


^crust 


dr 


The shear modulus in the crust, ignoring magnetic effects, is (Strohnrayer et al. 19911 
0.1194 m(Ze) 2 


(Bl) 

(B2) 

(B3) 

(B4) 


^ 1 + 0.595(173/T) 2 a ’ 

where rn is the number density of ions of charge Ze, a is the Wigner-Seitz cell radius given by rii4na 3 /3 = 1, and 
T = (Ze) 2 / (akT) where k is Boltzmann’s constant. Typically in the crust, T >> 173 and the second term in the denominator 


is negligible. For the composition of the inner crust, we use the results of Douchin & Haensel (2001), conveniently expressed 
analytically by |Haensel fe Potekhin| ( |2004[ ) . We solve for crust structure using the Newtonian equation for hydrostatic equi¬ 
librium, for a stellar radius of 10 km and a stellar mass of 1.4 Mq. In the evaluation of the shear-wave speed in the crust, we 
include the effects of nuclear entrainment (Charnel 20051 20121. Further details are given in Link (20141. 


We find that /i crua t drops from 2 x 10 3U erg cm” 3 at the base of the crust to 8 x 10 28 erg cm -3 over 370 m. Because 
we have assumed a homogeneous crust, we take the geometric mean of /i crU st over this range, which is 4 x 10 29 erg cm -3 , 
and set the crust thickness to 370 m. (The matter in the remainder of the crust contributes less to the rigidity than does the 
tangled field). The shear speed in the crust varies by a factor of about two in this region, with a geometric mean of 5 x 10 _3 c. 
With these values, we solve the boundary conditions numerically for the corrected eigenfrequencies. Crust rigidity increases 
the mode frequency by < 12% for each fundamental of a given l, and < 3% for harmonics. 

The treatment by Douchin & H aensel| ( [20 01) of the inner crust gives somewhat higher values of the shear modulus at 
the base of the crust than do other studies. The equation of state of . 


Aknral et al. 


the base of the crust that is about 0.6 the shear speed of Douchin & Haensel (2001), and a corresponding shear modulus that 


(1998), for example, gives a shear speed at 


is smaller by a factor of about 2.8. Reducing the crust shear modulus by a factor of two, and repeating the above analysis, 
shows that crust rigidity increases the fundamental frequency for each l by less than 10%. 
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